Invited talk at the "GAMMA2001" Symposium, April 4-6, 2001, Baltimore. 

To be published in "Gamma-Ray Astronomy 2001" (American Institute of Physics) 



Gamma-Ray Signatures of Supernovae and 

Hypernovae 

Ken'ichi Nomoto*, Keiichi Maeda*, Yuko Mochizuki^ Shiomi Kumagai**, 
Hideyuki Umeda*, Takayoshi Nakamura* and Isao Tanihata^ 

*Department of Astronomy and Research Center for the Early Universe, School of Science, 
University of Tokyo, Bunkyo-ku, Tokyo 113-0033, JAPAN 
^ RIKEN (The Institute of Physical and Chemical Research), Hirosawa 2-1, Wako, 

Saitama 351-0198, JAPAN 
**Department of Physics, College of Science and Technology, Nihon University, 
Kanda-Surugadai 1-8, Chiyoda-ku, Tokyo 101, JAPAN 

Abstract. 

We review the characteristics of nucleosynthesis and radioactivities in 'Hypernovae', i.e., su- 
pernovae with very large explosion energies (> 10^^ ergs) and their y-ray line signatures. We also 

discuss the ^"^Ti line y-rays from SN1987A and the detectability with INTEGRAL. Signatures of hy- 
pernova nucleosynthesis are seen in the large [(Ti, Zn)/Fe] ratios in very metal poor stars. Radioac- 
tivities in hypernovae compared to those of ordinary core-collapse supernovae show the following 
characteristics: 1) The complete Si burning region is more extended, so that the ejected mass of 
^''Ni can be much larger. 2) Si-burning takes place in higher entropy and more a-rich environment. 
Thus the ^^Ti abundance relative to ^^Ni is much larger. In aspherical explosions, '*'*Ti is even more 
abundant and ejected with velocities as high as ^ 15,000 km s^\ which could be observed in y-ray 
line profiles. 3) The abundance of ^^Al is not so sensitive to the explosion energy, while the ^''Fe 
abundance is enhanced by a factor of ~ 3. 



INTRODUCTION 

Massive stars in the range of 8 to ~ lOOM© undergo core-collapse at the end of their 
evolution and become Type II and Ib/c supernovae (SNe II and SNe Ib/c). These SNe II 
and SNe Ib/c release large explosion energies and eject explosive nucleosynthesis prod- 
ucts, thus being major sources of radioactive species. Until recently, we have considered 
supernovae with the explosion energies of £ = 1 - 1.5 x 10^^ ergs. These energies have 
been estimated from the observations of nearby supernovae, such as SNe 1987 A, 1993 J, 
and 19941. Also the progenitors of these SNe are estimated to be 13 - 20 Mq stars (e.g., 

[0)- 

Recently, there have been a number of candidates for the gamma-ray burst 
(GRB)/supernova (SN) connection (see ^ for references), including GRB980425/SNIc 
1998bw, GRB971115/SNIC 1997ef, GRB970514/SNIIn 1997cy, GRB980910/SNIIn 
1999E, GRB980326, and GRB970228. Among the SNe with a possible GRB counter- 
part, SNeIc 1998bw ^ and 1997ef ^ are characterized by a very large kinetic 



explosion energy, E > 10^^ erg. This is more than one order of magnitude larger than 
in typical SNe, so that these objects may be called "Hypemovae". These SNe produced 
more ^^Ni than the average core collapse SN. The masses of these hypernova progeni- 
tors are estimated to be M ^ 25Mq. These massive stars are likely to form black holes, 
while less massive stars form neutron stars (see, however, 

Regarding y-ray signatures of such hypernovae, whether and how the hypernovae 
actually induce gamma-ray bursts needs further study of aspherical explosions (e.g., 
[0, Hi). Another y-ray signatures, we discuss here, are the line y-ray emissions. We review 
the characteristics of nucleosynthesis and radioactivity in hypernovae and their y-ray line 
signatures. We also discuss the ^^Ti line y-rays from SN1987A and the detectability with 
INTEGRAL. For line y-rays from Type la supernovae, see Kumagai & Nomoto [Q] for 
a review. 

Before discussing y-rays, we first point out that signatures of hypernova nucleosyn- 
thesis are seen in the large [(Ti, Zn)/Fe] ratios in very metal poor stars. 



NUCLEOSYNTHESIS IN HYPERNOVA EXPLOSIONS 

In core-collapse supernovae/hypernovae, stellar material undergoes shock heating and 
subsequent explosive nucleosynthesis. Iron-peak elements are produced in two distinct 
regions, which are characterized by the peak temperature, Tpeak, of the shocked material. 
For Tpeak > 5 X lO^K, material undergoes complete Si burning whose products include 
Co, Zn, V, and some Cr after radioactive decays. For 4 x lO^K < Tpeak < 5 x lO^K, 
incomplete Si burning takes place and its after decay products include Cr and Mn (e.g., 

m E3 0)- 

We note the following characteristics of nucleosynthesis with very large explosion 
energies [[T3[]: 

1) Both complete and incomplete Si-buming regions shift outward in mass compared 
with normal supernovae, so that the mass ratio between the complete and incomplete 
Si-buming regions is larger. As a result, higher energy explosions tend to produce larger 
[(Zn, Co)/Fe], small [(Mn, Cr)/Fe], and larger [Fe/0]. The elements synthesized in this 
region such as ^^Ni, ^^Cu, ^^Zn, and ^^Ge (which decay into ^^Co, ^^Co, ^^Cu, and ^'^Zn, 
respectively) are ejected more abundantly than in normal supernovae. 

2) In the complete Si-buming region of hypernovae, elements produced by a-rich 
freezeout are enhanced because nucleosynthesis proceeds at lower densities (i.e., higher 
entropy) and thus a larger amount of ^He is left. Hence, elements synthesized through 
capturing of a-particles, such as '^^Ti, ^^Cr, and ^"^Ge (decaying into ^^Ca, ^^Ti, and 
^^Zn, respectively) are more abundant. 

3) Oxygen buming takes place in more extended, lower density regions for the larger 
explosion energy. Hence, O, C, Al are bumed more efficiently and their abundances in 
the ejecta are smaller, while a larger amount of buming products such as Si, S, and Ar are 
synthesized. Therefore, hypemova nucleosynthesis is characterized by large abundance 
ratios of [Si/0], [S/0], [Ti/0], and [Ca/0]. 
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FIGURE 1. The isotopic composition of the ejecta in the direction of the jet (upper panel) and 
perpendicular to it (lower panel). The ordinate indicates the initial spherical Lagrangian coordinate (Mr) 
of the test particles (lower scale), and the final expansion velocities (V) of those particles (upper scale) 



ASPHERICAL EXPLOSIONS 

Nakamura et al. [[T^ and Mazzali et al. [ p3| ] have identified some signatures of asym- 
metic explosion in the late light curve and spectra of SN 1998bw. Maeda et al. [ [T^ ] have 
examined the effect of aspherical (jet-like) explosions on nucleosynthesis in hypernovae. 
The progenitor model is the 16 Mq He core of the 40 Mq star and the explosion energy 
is £■ = 1 X 10^^ ergs. 

Figure |T] shows the isotopic composition of the ejecta of asymmetric explosion model 
in the direction of the jet (upper panel) and perpendicular to it (lower panel). 

In the z-direction, where the ejecta carry more kinetic energy, the shock is stronger 
and post-shock temperatures are higher. Therefore, larger amounts of a-rich freeze-out 
elements, such as ^He, ^^Ti, and ^^Ni are produced in the z-direction than in the r- 



direction. 

On the other hand, along the r-direction ^^Ni is produced only in the deepest layers, 
and the elements ejected in this direction are mostly the products of hydrostatic nuclear 
burning stages (O) with some explosive oxygen-burning products (Si, S, etc). 

In the spherical case, Zn is produced only in the deepest layer, while in the aspherical 
model, the complete silicon burning region is elongated to the z (jet) direction, so that 
[Zn/Fe] is enhanced irrespective of the mass cut. On the other hand, ^^Mn, which is 
produced by incomplete silicon burning, surrounds ^^Fe and located preferentially in 
the r-direction. 

In this way, larger asphericity in the explosion leads to larger [Zn/Fe] and [Co/Fe], 
but to smaller [Mn/Fe] and [Cr/Fe]. Then, if the degree of the asphericity tends to be 
larger for lower [Fe/H], the trends of [(Zn, Co, Mn, Cr)/Fe] follow the ones observed in 
metal-poor stars, as discussed later. 



SIGNATURES OF HYPERNOVA NUCLEOSYNTHESIS IN 
GALACTIC CHEMICAL EVOLUTION 

Several observational signatures of hypemova nucleosynthesis have been noticed in 



several objects [|13|]. The abundance pattern of metal-poor stars with [Fe/H] < —2 
provides us with very important information on the formation, evolution, and explosions 
of massive stars in the early evolution of the galaxy. 

In the early galactic epoch when the galaxy is not yet chemically well-mixed, [Fe/H] 
may well be determined by mostly a single SN event [|T^. The formation of metal-poor 
stars is supposed to be driven by a supernova shock, so that [Fe/H] is determined by the 
ejected Fe mass and the amount of circumstellar hydrogen swept-up by the shock wave 
[[T9|]. Then, hypemovae with larger E are likely to induce the formation of stars with 
smaller [Fe/H], because the mass of interstellar hydrogen swept up by a hypemova is 



roughly proportional to E [ [191 , |2^] and the ratio of the ejected iron mass to E is smaller 
for hypemovae than for canonical supemovae. 

The observed abundances of metal-poor halo stars show quite interesting pattem. 
There are significant differences between the abundance patterns in the iron-peak el- 
ements below and above [Fe/H]~ —2.5 — 3, which cannot be explained with the con- 
ventional chemical evolution model that uses previous nucleosynthesis yields. 

1) For [Fe/H]^ —2.5, the mean values of [Cr/Fe] and [Mn/Fe] decrease toward 
smaller metalicity, while [Co/Fe] increases [p1,[T^]. 



2) [Zn/Fe]~ for [Fe/H] ~ -3 to [El, while at [Fe/H] < -3.3, [Zn/Fe] increases 



toward smaller metalicity [|23|, |24|]. 

The larger [(Zn, Co)/Fe] and smaller [(Mn, Cr)/Fe] in the supemova ejecta can be 
realized if the mass ratio between the complete Si buming region and the incomplete 
Si burning region is larger, or equivalently if deep material from complete Si-burning 
region is ejected by mixing or aspherical effects. This can be realized if (1) the mass 



cut between the ejecta and the collapsed star is located at smaller M,. [|2q], (2) E is 



larger to move the outer edge of the complete Si burning region to larger My [|27| ], or (3) 
asphericity in the explosion is larger. 





FIGURE 2. Left: The distribution of Ni (open circles) and O (dots) with density contours (Hnes) 
[]r^. Right: The density distributions of ^^Ti (left half) and ^^Ni (right half). The densities of each element 
are represented with linear scale, from to the max density of each element pS]]. 



Also a large explosion energy E results in the enhancement of the local mass fractions 
of Zn and Co, while Cr and Mn are not enhanced [p5|]. Models with £51 = 10^^ ergs 
do not produce sufficiently large [Zn/Fe]. To be compatible with the observations of 
[Zn/Fe] ~ 0.5, the explosion energy must be much larger, i.e., > 20 for M > IOMq. 

Therefore, if hypernovae made significant contributions to the early Galactic chemical 
evolution, it could explain the large Zn and Co abundances and the small Mn and Cr 
abundances observed in very metal-poor stars. 



NUCLEAR RADIOACTIVITY IN HYPERNOVAE 



Hypernovae produce radioactive species in larger amount than normal SNe II, thus being 
important sources of line y-rays and positrons. 



^^Ni and Positrons 

A large amount of ^^Ni can be produced as actually observed: ~ 0.5 Mq from 
SN1998bw, ~ O.15M0 from SN1997ef, and even much larger from SN1999as. Thus 
hypernovae are the important source of positrons as well. In particular, positrons in the 
galactic center region might be significantly contributed by hypernovae, because as ob- 
served in M82 [p3|], hypernovae would make an important contribution in nucleosynthe- 
sis associated with starburst events. 



TABLE 1. The mass ratio ^^Ti/S^Ni for 40Mq models [||] . Values of 
spherical models are taken from Nakamura et al. [^7|]. 





Spherical (E^i ~ I) 


Spherical (E51 = 10) 


Aspherical (£51 = 10) 


6.51E-4 


2.79E-3 


l.OlE-2 



TABLE 2. Mass of ^^-Al and "'Fe (M0) (Z=0.02) 
& 





£51 


1 


10 


1 


10 


M 


26 


Al 


60 


Fe 


25Mq 
40Mq 


8.87e-5 
1.17e-4 
8.61e-5 


5.73e-5 
1.22e-4 
1.20e-4 


1.37E-5 
1.86E-4 
1.82E-5 


1.31E-4 
3.38E-4 
2.47E-5 



44Ti 

'^'^Ti is synthesized in the a-rich freezeout, in complete Si burning. In the very high 
entropy environment in hypernovae, the photons dissociate all the preexisting nuclei 
down to essentially a-particles and neutrons. The abundances then shift to those in 
nuclear statistical equilibrium but freezes out with excess a-particles. The ^^Ti yield 
depends strongly on the location of the mass cut, electron fraction, and entropy condition 
in the a-rich freezeout. Therefore the initial yield of ^^Ti provides unique information 
to constrain the explosion models. 

Iron-peak elements are produced in the deep region near the mass cut, so that their pro- 
duction is strongly affected by asphericity of an explosion. Figure |^ shows the 2D den- 
sity distribution of ^^Ti and ^^Ni, which is distributed preferentially in the z-direction. 
^'^Ti, which is produced by the strong a-rich freezeout, is distributed preferentially in the 
z-direction. Moreover, ^^Ti production is strongly enhanced compared with a spherical 
model since the post-shock temperature along the z-direction is much higher 



than that of a spherical model. 

Table |T] summarize the mass ratio ^^Ti/^^Ni produced in the spherical and the as- 
pherical explosion of the 40 Mq models. We can see that ^^Ti is strongly enhanced in 
hypernova and aspherical models The y-rays from the decays of ^^Ti are possible 
tools to investigate the Galactic hypernova remnants. 



2^A1 and ^"Fe 

Preliminary results of ^^Al and ^^Fe are summarized in Table ^ [^5]]. 

1) The ^^Al abundance does not depend much on £. It is a little smaller in hypernovae 
than supemovae, because more ^^Al is consumed in oxygen burning as are O, Mg, Ne. 

2) ^^Fe abundance is larger by a factor of ~ 1.5 - 10. 



RX J0852-4622/GRO J0854-4622 



COMPTEL has detected the ^^Ti 1157 keV y-ray line from the supernova remnant 
(SNR) RX J0852-4622 [PO, M, P], though the evidence is found at the 2o to 4a 



significance level [^3]] . ASCA has observed RX J0852-4622 and detected the 4. 1 keV X- 
ray emission line from Ca [|34l]. Tsunemi et al. [Q] have provided the following analysis 
and interpretation. The abundance of Ca is oversolar by a factor of 8 ± 5, while other 
elements are subsolar. The mass of Ca is ~ 1.1 X 10^3^0. The excess Ca is likely to be 
^^Ca, the decay product of "^"^Ti. This feature is seen only in the north-west shell, which 
suggests that the supernova ejecta has just collided with the interstellar material there. 

lyudin & Aschenbach p5| ] have assumed that the width of the ^^Ti line is due to 
Doppler broadening and thus "^"^Ti is expanding at ~ 15,000 km s^^ [Q. Then they 
suggested that such high velocity ^"^Ti is ejected from a Type la supernova of sub- 
Chandrasekhar mass because the model produces ~ 10^^ Mq ^Ti in the outer 

He detonation zone which expands at ~ 15,000 km s^^ 

Here we suggest an alternative model for the high velocity ^^Ti. As seen in Figure |^, 
the asymmetric hypemova explosion ejects ^^Ti at ~ 15,000 km s^^ in the jet direction. 
The amount of ^^Ti is as large as 1 x 10^-^Mq (Table |I]), being consistent with the 
observation. 

INTEGRAL observations of the "^"^Ti lines and their line profiles from RX J0852-4622 
are important to discriminates the models and clarify the energetics of the explosion. 



SN 1987A 

SN 1987 A in the LMC has shown for the first time that the energy source of supernova 
ejecta directly comes from the decays of radioactive nuclei (e.g., pB|, for reviews). 



In this section, we investigate the initial abundance of ^^Ti to discuss the detection 
possibility of the line gamma-rays from SN 1987A, by comparing theoretical light 
curves with the observed bolometric luminosity of SN 1987 A. 

In SN 1987 A, it is established that the observed light curve in early time is first 
governed by ^^Ni [ti/2 (half-life) = 6.1 d] and then its daughter ^^Co (ti/2 = 77.3 d). 

The synthesized ^^Co nuclide decays to stable ^^Fe. As we shall see later, the observed 
light curve in the wavelength ranging from ultraviolet (UV) to infrared (IR) has been 
successfully modeled with the energy supply from the decay of ^^Co until ~ 800 days 
(e.g., [^8|]). This has been directly confirmed by the detection of the hard X-rays and the 
line y-rays from the decay sequence of ^^Co. 

Afterwards, the decline of the observed light curve apparently slowed down, due to 
the decay of ^^Co {ti/2 = 272 d). The slowness of the decline of observed light curve 
becomes distinguished in particular after ~ 1500 days from the explosion. The dominant 
energy source at this moment can be attributed to ^"^Ti decay. ^^Ti decays by orbital 
electron capture to ^'^Sc, emitting 67.9 keV (100 %) and 78.4 keV (98 %) lines. "^^Sc 
then decays mainly by positron emission into "^^Ca, which emits a 1157 keV (100 %) 
de-excitation line. 

Obviously, a luminosity observation in SN 1987 A of late years is crucial to study the 



property of the extra energy source, namely, the ^"^Ti production. Recently, this crucial 
luminosity at 3600 days, i.e., 10 years after the explosion was reported by Suntzeff [HO]. 
The observed luminosity in the UV-IR range is L = (1.9 ±0.6) x 10-^^ergsec^^. The 
result is obtained by the collaboration of CTIO with HST, and hereafter we refer to this 
as the CTIO+HST luminosity. In the following, we investigate whether the ^"^Ti decay 
provides enough energy to account for the observed CTIO+HST luminosity under the 
latest half-life value of ^^Ti [P]. 



The Half-Life of 44Ti 

The energy release from the ^'^Ti decay depends strongly on its half-life. With recent 
experimental efforts, the half-life appears to be settled in 60 ± 3 years, including the 
errors up to 3 o (see [^] and references therein). 

We note here that the half-life values obtained in laboratories are for neutral atoms. 
Since ^^Ti undergoes pure orbital electron capture decay, its decay rate becomes smaller 
than the experimental value if '^'^Ti is highly ionized under the condition of a supernova. 
For example, the decay rates of hydrogen-like and helium-like ions are, respectively, ~ 
44 % and ~ 88 % of those of neutral atoms. Details are found in Mochizuki et al. [ ]43| ] 
and Mochizuki Q- 

It is expected that "^^Ti was neutral when the CTIO+HST observation was carried 
out. We thus adopt 60 ± 3 years to calculate theoretical light curves to compare the 
observation. 



Theoretical Light Curve of SN 1987A 

Our calculation of the light curves is based on Kumagai et al. [FSp and the adopted 



nuclear decay property data are updated. We perform Monte Carlo simulations of the 
Compton degradation of the line y-rays emitted from the decays of ^^Co (14 keV, 122 
keV, 136 keV, etc.) and ^^Ti (68 keV, 78 keV, 1157 keV) to obtain the UV-IR light 
curves. The UV, optical, and IR photons originate from the energy loss of the emitted 
y-rays during the radiative transfer in the ejecta. The calculated UV-IR luminosity is the 
result of subtracting the energy of the X-ray and y-ray photons which have managed to 
get out of the remnant. 

For the velocity distribution of particles, we adopt the explosion model 14E1 proposed 
by Shigeyama & Nomoto whose main-sequence mass, ejecta mass, and explosion 
energy are 20 Mq, 14.6 Mq (4.4 Mq core material plus 10.2 Mq hydrogen-rich enve- 
lope), and 1 X 10^^ erg, respectively. This model was derived from a detailed analysis of 
the plateau shape of the light curve of SN 1987 A which is observed until 120 days after 
the explosion, and well accounts for the earlier optical. X-ray, and y-ray light curves of 
SN 1987A g^. Note that the ^^Ni mass in SN 1987A has been constrained as 0.07 Mq 



from the intensity during the observed exponential decline. The synthesized mass of the 
radioactive ^^Ni is adopted from Hashimoto, Nomoto, & Shigeyama [ [lO|l , and the input 
value of the initial ^'^Ti abundance is varied within the uncertainty studied in Hashimoto 
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FIGURE 3. The calculated light curve and the observed bolometric (UV to IR) luminosity of SN 1987A, 
including the latest observed luminosity (CTIO+HST). 



et al. [0 to compare the CTIO+HST luminosity. 

In Figure ^ we show the calculated UV-IR light curve (solid) and the observed lu- 
minosities of SN1987A (CTIO and ESO). For this, the ^^Ti half-life of 60 yrs and 
<^^Ti/^^Ni>= 1 have been adopted. Here, <^^Ti/^^Ni> is defined as the ratio of 
^^Ti/^^Ni in SN1987A to ^^Ca/^^Fe in the solar neighborhood, i.e., <44Ti/56Ni> = 
[Z(44ji)/X(5^Ni)]/[X(44Ca)/X(5^Fe)]0. The three decay sequences of ^^Ni, ^^Ni and 

Ti are used for the energy sources as shown in Figure ^ We found that the observed 
CTIO+HST luminosity is reasonably explained with the energy release from the "^^Ti 
decay for its half-life of 60 ± 3 yrs. 



Production of 44Ti in SN 1987A 

The <44xi/56Ni> values that are required to explain the CTIO+HST luminosity are 
found to be roughly between 1 and 2. Since the synthesized mass of ^^Ni has been 
determined to be O.OVM© for SN 1987A, the obtained <'^'^Ti/^''Ni> values can directly 
be translated to the initial mass of ^"^Ti. Note that the amount of "^^Ti is the total value 
that is responsible for the overall energy in y-rays, X-rays, and in the UV-IR range. 

We thus obtained that the initial ^'^Ti mass is (1.1 - 2.5) X 10^4 foj- the half- 
life value of 60 ± 3 yrs. We remark that the ^^Ti yield calculated from recent explosive 
nucleosynthesis models overlap with our estimate (e.g., [ ^7] , 0). However, it should 
be also mentioned that these model predictions are subject to nuclear reaction cross 
sections which have not yet measured so far (see, [ |50| , |5T1]). 

Note that Chugai et al. [151] estimated that positrons from (1-2) x lO^'^M© of "^"^Ti 



provides the overall luminosity of the Fell emission lines, and Lundqvist et al. [53, 54] 
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FIGURE 4. Predictions of the line gamma-ray fluxes from the ^^Ti decay sequence (filled circles) for 
file initial mass of 44Ti of 1 x 10-'*MQand2x IO^^^Mq, respectively. The sensitivity limit of INTEGRAL 
SPI (10^ sec) is also shown by the dash-dotted curve. 



obtained the upper limit of the ^^Ti mass, 1 - 1.5 x 10"^ M©, based on ISO SWS/LWS 
observations. 
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Ti in SN 1987A with INTEGRAL 



Detection of ^^Ti line y-rays from SN 1987 A should be an important target for the 
EsTTEGRAL mission. In Figure ^ we show the expected fluxes from SN 1987A on 
Earth for 68, 78, and 1157 keV lines, respectively. For each y-ray line in Figure ^, the 
upper and lower filled circles are calculated for the '^'^Ti yield of 2 x 10^^ Mq, and 1 
X 10""^ Mq, respectively. In Figure the sensitivity limit of INTEGRAL SPI (10^ sec) 
is also shown, which strongly depends on energy. We see from Figure | that it is possible 
to marginally detect the 1 157 keV line if the produced amount of "^^i is as large as 2 
X 10^"* Mq, but for other cases the expected fluxes lie below the sensitivity limit of SPI. 

In the above discussion of the expected fluxes, no ionization effect on the ^^Ti decay 
is taken into account. Recently, helium-like and hydrogen-like ions of O, Ne, Mg, and 
Si have been observed with Chandra X-ray observatory [pS]]. The ionization of the 
ejected material was caused by the shock heating associated with the collision with 
the circumstellar matter (including the ring). As discussed previously, the decay rate of 
highly ionized ^^Ti can become considerably small compared with that of neutral "^^Ti. 
This means, as claimed by Mochizuki [p4|], all the ^^Ti y-ray line fluxes from SN 1987 A 
possibly lie below the detection limit of INTEGRAL SPI. 

Finally, we should note that INTEGRAL has capabilities both in the detection of 
nuclear lines and pulsed emission. It has not been clear yet whether SN1987A has 



formed a neutron star or a black hole. It is certainly worth trying fo INTEGRAL and 
Astro E-II to search for pulsation. 
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